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Abstract - We study graphene on a photonic crystal operating in the terahertz (THz) spectral 
range. We show that the absorption of graphene becomes a modulated function of frequency 
and can be enhanced by more than three times at specific frequency values, depending on the 
parameters of the system. The problem of a semi-infinite photonic crystal is also solved. 



The optical properties of graphene have been exten- 
sively studied both theoretically [I 14 and experimen- 
tally [15]j22] . More recently, the interest has shifted to 
how graphene interacts with electromagnetic radiation in 
the THz spectral range 23-28 . In particular, one of 



the goals is to enhance the absorption of graphene. This 
maybe done in different ways: (i) producing micro-disks 
of graphene on a layered structure [23] ; (ii) exploiting the 
physics of quantum dots on graphene [25] ; (iii) using a 
graphene based grating 29 30 ; (iv) putting graphene in- 
side an optical cavity [31 



32 



; and (v) depositing graphene 
on a photonic crystal, operating in the visible range of the 
spectrum UA]. In cases (i), (ii), and (iii) the excitation of 
plasmons is responsible for the enhancement of the absorp- 
tion. In case (iv) the radiation performs many round trips 
inside the cavity enhancing the chances of being absorbed 
by graphene. In case (v) the authors use a photonic crystal 
made of SiC^/Si. In the visible range of the spectrum the 
dielectric constants of Si02 and Si differ by more than one 
order of magnitude and choosing the width of the SiC^/Si 
appropriately it is possible to induce a large photonic band 
gap in the visible range. Combining the presence of the 
band gap with an initial spacer layer the absorption can 
be enhanced by a factor of four. In the case studied in 
Ref. 1 14 the optical conductivity of graphene is controlled 
by inter-band transitions. In this regime, the conductivity 
is frequency independent and reads 15 



2h 
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a value dubbed the universal conductivity of graphene. 

In general, the conductivity of graphene is a sum of 
two contributions: (i) a Drude term, describing intra-band 



processes and (ii) a term describing inter-band transitions. 
At zero temperature the optical conductivity has a simple 
analytical expression [1,2,5,9,10 . For what concerns our 
study, the physics of the system is dominated by the Drude 
contribution, which reads 
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where Y is the relaxation rate, ep > is the Fermi level 
position with respect to the Dirac point, and uo is the fre- 
quency of the incoming radiation. It should be noted that 
gb has a strong frequency dependence and is responsible 
for the optical behaviour of graphene in the THz spectral 
range. 

In what follows, the absorption of graphene on a system 
such as the one represented in Fig. [I] is computed. The 
micro-structure has a first layer (spacer) of Si02, of length 
d s , followed by a periodic array of Si/Si02 which extends 
for several layers. We have found that no more than few 
layers are necessary for obtaining an enhancement of the 
absorption. We write d s as 



d. 
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where j3 is a numerical parameter chosen to optimize 
the absorption, e s is the dielectric constant of the spacer 
(SiC^), and A is a parameter of the order of few microns 
also chosen to optimize the absorption. The lengths dj of 
the Si/Si02 periodic structure layers are given by 



dj 



A 
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where j = a, 6, with e a (ei ) ) standing for the dielectric con- 
stant of Si(Si02). In the THz we approximate the di- 
electric constants of Si and Si02 by their static values, 
e a = 11.9 and = 3.9. 




Fig. 1: Graphene on a micro-structure made of SiCb/Si. The 
initial spacer has a different length, d s , from the other layers 
of Si02 and Si, which have length dj. Graphene can be doped 
either using a transparent top gate or chemically [33 34 ; in the 
latter case no gate is necessary. 

To determine the absorption we have to compute the 
amount of light reflected and transmitted by and through 
the structure, respectively. If we denote the transmittance 
by T and the reflectance by 7£, the absorbance is defined, 
assuming normal incidence, as 



A=1-T-1Z. 
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It is also assumed that the radiation impinges from the air 
and emerges the structure to the air as well. Since we are 
considering normal incidence, no distinction is necessary 
between s— and p— polarized waves. 

An electromagnetic wave impinges on graphene from 
the left along the x— direction. Due to reflection on the 
structure the electric field has the form 

E = (E 0+ e ikx + E -e- ikx )u v . (6) 

From Maxwell's equations, the magnetic field is given by 



iujB z = d x E v 



(?) 



At the graphene interface the boundary conditions are 
El = Er and Bl — Br = hq(Jl>El, where the subscript 
L(R) stands for the fields to the left (right) of graphene. At 
the other interfaces the boundary condition for the mag- 
netic field is simply Bl = Br. On the right hand side of 
the structure we write the electric field as 



E 



771 ikx.. 
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Using the transfer matrix method we can relate the fields 
on the left of the structure to the fields on the right as 



, (9) 
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from where it follows 
T = 
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and 
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For the structure we are considering, the transfer matrix 
T is given by 



f=(K 0tb + afJ)^ 8 K bia t 



N 
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where N is the number of times the matrix T is repeated 
in the structure and T reads 



T = ^ a K a ^ b K t 



(13) 



The several matrices appearing in T and f are given 
below. The matrix characterizing the free propagation 
through a given material is 



e -iuj^ejdj/c 







Juy/ejdj/c 
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with j = s, a, b. The matrix referring to the transfer 
through the boundary of two dielectrics is given by 



Ki. 



(VV^+I)/ 2 (1-v^M)/ 2 

(l-V^)/2 (^M + l)/2 



(15) 



In particular, K a ^ and Kq^ refers to the passage from the 
air to the last slab of Si and from the first Si02 spacer to 
the air, respectively [thus eo = 1 in Eq. (15)]. The matrix 
Ko,b + &fJ characterizes the transfer through graphene to 
air. Here matrix 
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a « 1/137 is the fine structure constant and / is a dimen- 
sionless function defined by 



/ 
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We note that in the infinite structure we can apply Bloch 
theorem, and relate the electric field amplitudes in the 
unit cell m with those in the unit cell m + 1 using the 
transfer matrix T as 
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E r > 
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where q is the Bloch wave vector. Then the spectrum of 
the infinite structure can be obtained from the condition 
det|J — Te iq ^ da+dh ^ \ = (here I is the unit matrix). Using 
the property of transfer matrix det|T| = 1, we obtain that 
the spectrum of the infinite photonic crystal, 



Tr(f) = 2cos[(d + d 6 )(7], 



(19) 



is defined by trace of T only. From Eq. (19) the pho- 
tonic band gaps can be obtained with no difficulty. We 
can also access the photonic band gaps by transmittance 
calculations, a route we choose below. 
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Fig. 2: Intensity plot of the absorbance. Top panel: as func- 
tion of the parameter f3 and the photon energy ftw (in meV). 
Bottom panel: as function of A (in /mi) and huo. High inten- 
sity spots corresponds points in the parameter space where the 
absorbance attains its higher values. We have used N = 4, 
f3 = 0.85, A =20 /mi, HF = 2.5 meV, and e F = 0.25 eV, a 
typical value. For reference, 4.1 meV corresponds to 1 THz. 



In Fig. [2] we represent the absorbance as function of: (i) 
the parameter j3 and the photon energy two (in meV)[top 
panel]; (ii) the parameter A (in jim) and hu [bottom 
panel] . Recall that the parameter f3 controls the length of 
the Si02 spacer whereas A influences the length of the lay- 
ers constituting the periodic structure. The high-intensity 
spots correspond to regions in parameters space where the 
absorbance is higher. The axes in Fig. [2] intercept at one 
of those high intensity spots. The width of the high inten- 
sity regions along the hu direction gives the width of the 
resonance. The position of the high-intensity regions is 
sensitive to the parameter A, as can be seen in the bottom 
panel of Fig. [2] These type of plots allow the optimiza- 
tion of the structure in order to maximize the absorbance. 
An optimized absorbance peak is given in Fig. [3J The 



parameters chosen are those corresponding to the inter- 
ception of the horizontal and vertical lines in Fig. [2j In 
Fig. [3] we represent the absorbance, the reflectance, and 
the transmittance as function of the frequency of the in- 
coming photon. Since the dielectric constants we have 
used for Si02 and Si are real, all the absorption comes 
of graphene alone. The intricate behaviour seen on the 
reflectance and transmittance curves, on the other hand, 
is mainly due to the structure on the right hand side of 
graphene. The dashed lines in the same figure represent 
the behaviour of free standing graphene. 

Several aspects are worth mention. Although we have 
used a fairly small value of TV, a stop gap due to the pho- 
tonic band gap of the infinite crystal is already apparent. 
This happens because the difference between the dielec- 
tric constants of SiC>2 and Si is significant. Most impor- 
tant however is the enhancement of the absorbance near 
the center of the photonic band gap. We see that around 
an energy Huj « 15.3 meV the absorbance is almost four 
times larger than the absorbance of free standing graphene 
(assuming it can be doped some how, e. g., by chemical 
doping 33,34 ; if we had chosen graphene on a Si02 wafer 
of the size of the photonic crystal the absorbance would 
be smaller than that for the free standing case). If we 
had chosen somewhat different parameters the resonance 
would shift to the edges of the stop gap. In the absorbance 
panel we plot the cases N = 1, N = 4, and N = 8. As it 
can be seen, for N = 1 the absorbance is not as intense 
as in the other two cases. At the same time the main dif- 
ference is the appearance of new satellite peaks to the left 
and to the right of the main one when we change from 
TV = 1 to N = 4, and to N = 8. The nature of these 
peaks can be expressed in the following simplified man- 
ner. In the infinite periodic structure allowed bands are 
continuous functions of Bloch wave vector q, while in the 
periodic structure containing a finite number of periods 
N the effect of quantization takes place: allowed bands 
consist of the discrete set of the frequencies, which can 



be determined from Eq.(19) using the discrete Bloch vec- 
tors q n = n7r/(N + 1) for n G [1,7V] (see Ref. [35]). In 
particular, for N = 4 the mode with n = 4 correspond to 
frequencies uo ~ 11.4 meV and uj ~ 19.6 meV (here we take 
into account the frequency range of Fig j3] only), while for 
N = 8 the allowed frequencies are 11.2 meV (mode with 
n = 7), 12.3 meV, and 18.8 meV (mode with n = 8). In 
both cases the predicted frequencies correspond well with 
the absorption peaks below and above the stop gap seen 
in Fig(3j 

It is instructive to compute the absorbance, reflectance, 
and transmittance when the photonic crystal is semi- 
infinite. In this case the oscillations due to the finiteness 
of the crystal seen in Fig. [3] are absent and it is possible to 
derive simple analytical formulae for the above quantities. 



From Eq.(18) it follows that 
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Fig. 3: Absorbance (left), reflectance (center), and transmit- 
tance (right) as function of frequency. The dashed line repre- 
sents the same quantity but for free-standing graphene. We 
have used N = 8, ft = 0.85, A =20 /xm, hF 2.5 meV, and 
€f — 0.25 eV. The cases N — 1,4 are also shown in the ab- 
sorbance panel. 
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Fig. 4: Absorbance (left), reflectance (center), and transmit- 
tance (right) as function of frequency for a semi-infinite pho- 
tonic crystal. The dashed line represents the same quantity 
but for free-standing graphene. We have used /3 = 0.85, A =20 
/xm, HT — 2.5 meV, and ef — 0.25 eV. 



where are the matrix elements of the matrix T. Now 
the amplitudes for x < relate to the amplitudes in the 
first Si slab as 
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from which it follows that 



E 0+ = (M 11 +p(q)M 12 )El + , 
E - = (M 2 i + p{q)M 22 )E 1 a+ , 

leading to the reflectance amplitude 

Eq- M 2 i + p(q)M 22 
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from which the reflectance is computed as 1Z = |r| 2 . Be- 
cause in each slab there is always an incoming and an 
outgoing wave the calculation of the transmittance is less 
obvious. We have to compute two field ratios: 



and 
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E 0+ Mn + p(q)M 12 
Then, the transmittance is given by 
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(25) 



(26) 



(27) 



where the factor ^/e^ comes from the calculation of the 
flux. In Fig. [4] we plot the absorbance (left), reflectance 



(center), and transmittance (right) as function of fre- 
quency for a semi-infinite photonic crystal. The resonance 
in the absorbance coincides exactly with that seen in Fig. 
[3j but the oscillations also seen in Fig. [3] are not present 
here. 

In conclusion, it is possible to use a photonic crystal- 
like structure to enhance the absorption of radiation by 
graphene in the THz spectral range at specific frequen- 
cies. The effect depends on the geometric properties of 
the crystal. At the center of the photonic band gap the 
absorption increases at the expenses of the reflectance (the 
resonance can shift its position depending on the param- 
eters chosen). If the Fermi energy increases so does the 
absorption and the width of the peak. It is also clear that 
for a finite crystal the absorbance becomes a modulated 
function of frequency. We have chosen the widths of the 
Si02 and Si having the same dependence on A (exception 
made to the initial layer which has the extra parameter f3). 
If we had chosen them to be different we would have an ex- 
tra degree of freedom to tweak the absorbance curves, but 
the optimization procedure would become cumbersome. If 
we had sandwiched graphene between two dielectrics the 
absorbance would decrease because of the optical conduc- 
tivity would appear in the equations divided by the square 
root of the dielectric constant of the medium. For the pa- 
rameters of Fig. [3] the next stop gap appears at an energy 
of about 47 meV (not shown). There the enhancement of 
the absorbance is larger by a factor of four relatively to the 
free standing case. Finally we note that the device does 
not necessarily require a gate for doping graphene. Sta- 
ble hole doping by chemical methods have recently been 
achieved, without affecting the transparency of the mate- 
rial [33] . 
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